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Abstract. Given A £ 2, mxn and 6 £ Z m , we consider the issue of 
existence of a nonnegative integral solution x £ N" to the system of 
linear equations Ax — b. We provide a discrete and explicit analogue of 
the celebrated Farkas lemma for linear systems in R" and prove that 
checking existence of integral solutions reduces to solving an explicit 
linear programming problem of fixed dimension, known in advance. 



1 Introduction 

Let A £ Z mxn ,fe £ Z m and consider the problem of existence of a solution 
x £ N™ of the system of linear equations 

Ax = b, (1) 

that is, the existence of a nonnegative integral solution of the linear system 
Ax = b. 

Contribution. The celebrated Farkas Lemma in linear algebra states that 

{x £ R™ | Ax = b} ^ [u £ R m and A'u > 0] =>• b'u > (2) 

(where A' (resp. b') stands for the transpose of A (resp. b)). 

To the best of our knowledge, there is no explicit discrete analogue of 0. In- 
deed, the (test) Gomory and Chvatal functions used by Blair and Jeroslow in [5] 
(see also Schrijver in Corollary 23.4b]) are defined implicitly and recursively, 
and do not provide a test directly in terms of the data A, b. 

In this paper we provide a discrete and explicit analogue of Farkas Lemma for 
to have a solution x £ N™. Namely, when A and b have nonnegative entries, 
that is, when A £ N mxn ,b £ N™, we prove that Q has a solution x £ N" if 
and only if the polynomial z i— ► z b — 1 (:= z^ 1 ■ ■ ■ z h ^ — 1) of R[zi, . . . , z m ], can 
be written 



- 1 = 



3=1 



1) 



(3) 



for some polynomials {Qj} in R[zi, . . . , z m ] with nonnegative coefficients. In 
other words, 



{x £ N" | Ax = b} ^ <S> z & - 1 = ^Qj(z)(z A * - 1), (4) 

for some polynomials {Qj} in R[^i . . . . , z rn ] with nonnegative coefficients. (Of 
course, the if part of the equivalence in is the hard part of the proof.) 

m m 

Moreover, the degree of the Qj's is bounded by 6* := bj — min Ajk- 

3=1 j=i 

Therefore, checking the existence of a solution x £ N" to Aa; = b, reduces to 
checking whether or not there is a nonnegative solution y to a system of linear 
equations where (i) y is the vector of unknown nonnegative coefficients of the Qj 's 
and (ii), the (finitely many) linear equations identify coefficients of same power 
in both sides of ®. This is a linear programming (LP) problem with ns(b*) 
variables and s(b* + maxj . Ajk) constraints, where s(u) := ( m ^ u ) denotes 
the dimension of the vector space of polynomials of degree u in m variables. In 
addition, all the coefficients of the associated matrix of constraints are all or 
±1. For instance, checking the existence of a solution x £ N™ to the knapsack 
equation a'x = b, reduces to solving a LP problem with n(b + 1 — mffij- dj) 
variables and b + 1 + maxj a,j — mim, a,j equality constraints. This result is also 
extended to the case where A £ Z mxn , b £ Z m , that is, when A and b may have 
nonnegative entries. 

We call (0J a Farkas lemma because as ©, it states a condition on the dual 
variables z associated with the constraints Ax = b. In addition, let z := e A and 
notice that the basic ingredients b'X and A'X of (J5J, also appear in J3J via z b 
which becomes e b A and via z Aj which becomes e^ A . Moreover, if indeed z b — 1 
has the representation then whenever A £ R m and A'X > (letting z := e A ) 



e bX -l = ^Qj(e x \...,e x « 



e (A'\) 3 



> 



(because all the Qj have nonnegative coefficients), which implies b'X > 0. Hence, 
we retrieve that b'X > whenever A'X > 0, which is to be expected since of 
course, the existence of nonnegative integral solutions to {Q implies the existence 
of nonnegative real solutions. 

Methodology. We use counting techniques based on generating functions 
as described by Barvinok and Pommersheim in [5] and by Brion and Vergne in 
| lloj . to easily obtain a simple explicit expression of the generating function (or, 
Z-transform) F : C m -> C of the function / : Z m -> N, b ^ f(b), that counts 
the lattice points x £ N™ of the convex polytope fi := {x £ R" | Ax — b}. Then 
/ is the inverse Z-transform of F and can be calculated by a complex integral. 
Existence of a solution x £ N" to (Q is equivalent to showing that f(b) > 1, and 
by a detailed analysis of this complex integral, we prove that Q is a necessary 
and sufficient condition on b for /(£>) > 1. 



2 Notation and preliminary results 



For a vector b G R m and a matrix A G R mx ", denote by b' and A' G R™ xm 
their respective transpose. Denote by e m G R m the vector with all entries equal 
to 1. Let R[xi, . . . ,x n ] be the ring of real- valued polynomials in the variables 
x%, . . . , x n . A polynomial / G R[xi, . . . , x n ] is written 

x -» f( X ) = J2 f* xa = E f«xr 

a£N" a£N" 

for finitely many real coefficients {/ Q }. 

Given a matrix A G Z mxn , let A,- G Z m denote its j-th column (equivalently, 
the j-th row of A'); then for every z G C m , z" 4 ^ stands for 

Z A 3 ._ z ^i... z A mj _ e <Aj,ln 2 ) _ e (A'ln Z )^ 

If G N m then z Aj is a monomial of R[zi, . . . , z m \. 



2.1 Preliminary result 

Let A G Z mxrl , 6 G Z m and consider the system of linear equations 

Ax = b; x G N", (5) 

and its associated convex polyhedron 

Q := {x G R" | Ax = b; x > 0}. (6) 

It is assumed that the recession cone {x G R" | Ax — 0; x > 0} of i7, reduces to 
the singleton {0}, so that Q is compact (equivalently, fl is a convex polytope). 

By a specialized version of a Farkas Lemma due to Carver, (see e.g. Schrijver 
in [SJ (33), p. 95]), this in turn implies that 

{A G R m | A'X > 0} ^ 0. (7) 

Denote by b i— > f(b) the function / : Z m — > N that counts the nonnegative 
integral solutions x G N m of the system of linear equations JSJ, that is, the 
lattice points x G N" of 12. In view of Q, /(&) is finite for all b G Z m because 
fl is compact. Let F : C m — > C be the two-sided Z-transform of /, that is, 

z ~ F(z) := J2 f( u ) z ' U = E f{u)z^---z-^ (8) 

when the above series converges on some domain D C C m . It turns out that 
F(z) is well-defined on 

D:={zgC™| ^••■•4^'|>1 j = l,...,n}. (9) 



Proposition 1. Let A e Z mxn ,b E Z" and assume that ^Jjjj holds. The 

m 



for all z G Z m that satisfy 



|^| = \z^---zi^\ > 1 j = 1, 



(10) 



(11) 



Moreover, 
f(b) 



,b-e„ 



(2^)- y kl , =71 y |Zml=7mn ™ ( i_ z - 



■dz (12) 



■ dz, 



pm)* y zer n;=i(i-^- 

wra£/i f := {z 6 C m | |zj = 7j}, and where 7 6 R™ is /ized and satisfies 

l Aj = lt l3 ---l£r > 1 j = l,...,n. (13) 

Proof. The proof is a verbatim copy of that of Lasserre and Zeron in [7| where 
the linear system Ax < b (instead of Ax = b) was considered, but for the sake 
of completeness we reproduce it here. Apply the definition (JHJ of F to obtain : 



n*) = E 



ueZ" 



E 1 



xGN n , Ax—u 



E 

ieN" 



E ^ 



Now observe that 

-(Aar)i „-(Aa:) 2 



n i^ Aik -- 



k=l 



no 

fc=i 



Hence, when (II 11 1 holds we obtain 



n 00 



w = n e c« 



fc— 1 Xfe— 



fc=i 



which is l)l()[l. and (|12|l is obtained by a direct application of the inverse Z- 
transform (see e.g. Conway in [^j). It remains to show that, indeed, the domain 
defined in 111|) is not empty. But this follows from (JJJ. Indeed take z^ := e Afc for 
all = 1, . . . , m, for any A that satisfies J7J). 



3 Main result 

Before proceeding to the general case A s Z mx ", we first consider the case 
A G N mXTl where A (and thus 6) has only nonnegative entries. 



3.1 The case A E N™ x 



In this section A G N mx ™ and thus, necessarily b € N m (otherwise fl = 0). 

Theorem 1. Let A G N mx ",6 G N m . Then the following two statements (i) 
and (ii) are equivalent : 

(i) The linear system Ax — b has a solution x S N™. 

(ii) The real-valued polynomial z i— » z b — 1 := z bl • • • z b ™ — 1 can be written 

n 

z b -l = ^2Qj{z){z^-l) (14) 
j=i 

for some real-valued polynomials Qj G R[zi, . . . , z m ], j = 1, . . . ,n, all of them 
with nonnegative coefficients. 

In addition, the degree of the Qj 's in \1J$ is bounded by 

m m 

b* := » niin^.l,,. (15) 

J = l 3=1 

For a proof see 21 
3.2 Discussion 

(a) Let s(u) := ('"^") the dimension of the vector space of polynomials of degree 
u in m variables. In view of Theorem ^ and with 6* as in l|15|) . checking the 
existence of a solution x G N" to Ax — b reduces to checking whether or not 
there exists a nonnegative solution y to a system of linear equations with : 

— n x s(b*) variables, the nonnegative coefficients of the Qj's. 

n 

— s(b* + maxN Ajk) equations to identify the terms of same power in both 

k — 4 

3=1 

sides of ((H)) . 

This in turn reduces to solving a LP problem with ns(b*) variables and s(b* + 
maxfc J2j Ajk) equality constraints. Observe that in view of l|14|) . this LP has a 
matrix of constraints with only and ±1 coefficients. 

(b) In fact, from the proof of Theorem ^ it follows that one may even 
enforce the weights Qj in l|14|) to be polynomials in Z[z\, . . . , z m ] (instead of 
R[zi, . . . , z m \) with nonnegative coefficients (and even with coefficients in {0, 1}). 
However, (a) above shows that the strength of Theorem ^ is precisely to allow 
Qi G R[ z\, . . . , z m ] as it permits to check feasibility by solving a (continuous) 
linear program. Enforcing Qj G Z[zi, .. . , z rn ] would result in an integer program 
of size larger than that of the original problem. 

(c) Theorem Q] reduces the issue of existence of a solution x G N" to a 
particular ideal membership problem, that is, Ax = b has a solution x G N" if and 
only if the polynomial z b — 1 belongs to the binomial ideal I — (z Aj — l)j=i,.. n C 
R[zi, . . . , z m ] and for some weights Qj all with nonnegative coefficients. 



Interestingly, consider the ideal J C R[zi, . . . , z m , yi, . . . , y n ] generated by 
the binomials z Aj — yj, j = 1, . . . ,n, and let G be a Grobner basis of J. Using 
the algebraic approach described by Adams and Loustaunau in T| §2.8], it is 
known that Ax = b has a solution x £ N™ if and only if the monomial z b is 
reduced (with respect to G) to some monomial y a , in which case a £ N™ is a 
feasible solution. Observe that this is not a Farkas lemma as we do not know 
in advance a £ N™ (we look for it!) to test whether z b — y a £ J. One has to 
apply Buchberger's algorithm to (i) find a reduced Grobner basis G of J, and 
(ii) reduce z b with respect to G and check whether the final result is a monomial 
y a . Moreover, note that the latter approach uses polynomials mn + m (primal) 
variables y and (dual) variables z, in contrast with the (only) m dual variables 
z in Theorem ^ 



3.3 The general case 

In this section we consider the general case A £ Z TOX " so that A may have nega- 
tive entries. The above arguments cannot be repeated because of the occurence 
of negative powers. However, let a £ N™, £ N be such that 

Aj k ■= A jk + a k > 0; bj := bj + f3 > 0; k = 1, . . . ,n; j = 1, . . . , m. (16) 

Note that once a £ N n is fixed as in l)16[l. we can choose (3 £ N as large as 
desired. Moreover, as J? defined in JHJl is compact, we have 

n n 

max {S oijXj \Ax = b} < max{> ctjXj \ % £ ='■ P*( a ) < °°- (17) 

j=l 3=1 

Given a G N™, the scalar p*(a) is easily calculated by solving a LP problem. 
Next, choose p*(a) < (3 £ N, and let A £ N mxn ,b £ N m be as in flfy. Then 
the existence of solutions x £ N" to Ax — b is equivalent to the existence of 
solutions {x, u) £ N n x N to the system of linear equations 




b 

P- 



(18) 



3=1 

Indeed, if Ax = b with x £ N™, then 



Ax + e m 22 a 3 x j ~ e m a j x j = b + ([3 ~ f3) e r , 

3=1 3=1 



or, equivalently, 



3=1 



and thus, as f3 > p*(a) > 5Z^i a j x i ( c ^- we see * na ^ ( x i u ) with (3 



En 
.7=1 "j 



G N, is a solution of (fT%| . Conversely, let (x, u) £ N" x N be a 



solution of (|18|l . Then, using the definitions of A and b, 

n n 

Ax + e m ttjXj +ue m = b + (3 e m ; Qj'^j + w = A 

so that Ax = b. The system of linear equations l|18l) can be put in the form 

A I e r , 



D 



X 




T 


u 







with B := 



a' I 1 



(19) 



and as B G pj(m+i)x(n+i) 5 we are back to the case analyzed in A3. II 



Theorem 2. Let A G z mx ™, b G Z™ and assume that fl defined in ^Bjj is com- 
pact. Let A S N mx ", b 6 N m , a G N" and (3 G N oe as in with (3 > p*(a) 
(cf. ). Then the following two statements (i) and (ii) are equivalent : 

(i) The system of linear equations Ax = b has a solution x G N™ . 

(ii) The real-valued polynomial z t—t z b (zy)^ — 1 G R[zi, ■ • ■ , z mi y] can be 
written 



z\zyf-\ = Q Q (z,y)(zy-l)+J2Qj(^y)( zAi ^y) ai - 1 ) ( 20 ) 

j=i 

for some real-valued polynomials {<5j}" =0 ^ n R[ z ii • • • > z m, y]> oil with nonnega- 
tive coefficients. 

In addition, the degree of the Qj 's in \2(J\) is bounded by 



(m + \)[3 + bj — min 



m + 1 , min 

k— l,...,n 



(m + l)a k + ^2 Ajk 
j'=i 



Proof. Apply Theorem ^ to the equivalent form i|19fl of the system Q in l|18[l 
(since B G N(' Tl+1 ) x ( n+1 ) and (b, /?) G N m+1 ), and use the definition (O of 
(6, /3) and 2. 



4 Proof of Theorem [T] 

Proof, (ii) =>• (i) . Assume that z b — 1 can be written as in l|14|) for some polyno- 
mials {Qj} with nonnegative coefficients {Qj a }, that is, 

Qj( z ) = Qja^ = QjaZ^-'-Z^, 

a6N m a€N m 



for finitely many nonzero (and nonnegative) coefficients {Qja}- By Proposition 
^ the number /(&) of nonnegative integral solutions x £ N" to the equation 
Ax = b, is given by 

If f z b-e m 



(27ri) m y kl | =71 j\ Zml = lm iiuo-- z ~ J] 

Writing z b ~ Cm as z~ em (z b — 1 + f ) we obtain 

f(b) = 5i+B 2 , 

with 



and 



(2^)™ i,, ll=71 i| Zro | =7m nLi(i-^ 



V V Qj ° / ... / — : dz 



From l|12[l in Proposition ^ (with 6 := 0) we recognize in B\ the number of 
solutions x £ N™ to the linear system Ax = 0, so that B\ = 1. Next, again from 
l|12l) in Proposition ^ (now with b := Aj + a), each term 

Qja f f z A j +a~e m 



(2m)™ J M=71 J\ Zm]=lm n^(l - 
is equal to 

Q 3Q x the number of integral solutions x £ N™ _1 

of the linear system A^x = Aj + a, where A^> is the matrix in N mx (™-i) 
obtained from A by deleting its j-th column Aj. As by hypothesis, each Qj a is 
nonnegative, it follows that 

n 

B ^ = E E c ^ * °> 

j=l aEN" 1 

and so f(b) = B\ + B2 > 1. In other words, the sytem Ax = b has at least one 
solution x £ N™. 

(i) => (ii). Let x £ N" be a solution of Ax = b, and write 



1 = z AlX1 - 1 + z AlX1 {z A2X2 -!) + ■■■ + *5-£=i A ^ (z 4 " 1 " - 1), 



and 



z A jXj _ 1 = ^A, _ y 



1 + z Aj H h z A ^ X] 15 



j = l,...,n, 



to obtain l|14(l with 



Qj(z) := zSLi [1 + + • • • + z A ^ x i-^ 



J = 2, • 



and 



z i ► Qx(z) := 1 + z Al H 



„A!(xi-l) 



We immediately see that each has all its coefficients nonnegative (and even 
in {0,1}). 

Finally, the bound on the degree follows immediately from the expression of 
the Qj's in the proof of (i) => (ii). 
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